Assessment of groundwater quality for human consumption and its health risks in the Middle Magdalena Valley, Colombia

Groundwater is the primary source of water for people living in rural areas, especially during seasons when surface water is contaminated or unavailable. In Colombia, people use groundwater as drinking water without additional treatment. In addition, there is no infrastructure for wastewater collection and sewage treatment in the region of the Middle Magdalena Valley. The current study aims to evaluate the quality of groundwater in this region to determine any potential health hazards associated with its consumption. To reach the objective, three (3) physicochemical and microbiological sampling campaigns were carried out during different hydrological periods. A total of 428 groundwater samples were analyzed for over 28 parameters. The results were compared with the water quality standards proposed by the US EPA and Colombian regulations for human consumption. The analysis revealed the presence of total and fecal coliforms in 89% and 58% of the analyzed samples, respectively, identifying them as the main contaminants in groundwater. Furthermore, the pH levels did not meet the standards set by the US EPA in 33.8% of the cases and by Colombian regulations in 31.02%. Additionally, 32.8%, 17.6%, 14.3%, and 10.9% of the samples failed to meet the established thresholds for apparent color, magnesium, iron, and nitrates, respectively, under both standards. Moreover, only the analyses of selenium, mercury, and zinc complied with the quality standards under both regulatory frameworks. Based on the Colombian Drinking-Water Quality Risk Index (CDWQRI-IRCA), the risk associated with water quality meant for human consumption was assessed. The results showed that over 84% of the samples analyzed posed a high risk to human health, 4.6% posed a medium risk, 5.5% posed a low risk, and only 5.7% posed no risk at all. Additionally, official mortality statistics for children under four years old were reviewed, which revealed two deaths in 2019 due to Acute Diarrheal Disease (ADD) caused by consumption of contaminated water. Therefore, it is crucial to implement water treatment systems, establish aqueducts in rural areas, and conduct rigorous and systematic monitoring of drinking water to ensure it is safe for human consumption. It is also important to track morbidity and mortality rates associated with water consumption.


Study area
The MMV is located in the northeastern region of Colombia, a plain zone with hills that do not exceed 288 masl (Fig. 1a).It has a surface area of 8500 km 2 and is geographically situated between latitudes 8°30ʹ00″ N-6°54ʹ00″ N and longitudes 73°21ʹ00″ W-74°06ʹ00″ W. VMM is bounded by the piedmont of the Central and Eastern Ranges of the Andes (Fig. 1b).
The mean temperature is 27°.The rainfall regime has a bimodal behavior characterized by two rainy periods.The first one between March and May and, the second one between August and November.The annual average rainfall is about 2000 mm/y, showing a decreasing tendency from the southern (3700 mm/y) to the northern area (1000 mm/y) 40,44 (Fig. 1c).The MMV is highly affected by the El Niño-Southern Oscillation (ENSO).During the wet season, floods affect communities and large areas of crops for long periods; in contrast, it is estimated a decrease of 60% in the run-off in the Magdalena-Cauca basin on a typical dry year 40,45,46 .This situation would be worse due to the increase in the frequency of extreme events associated with climate change.Even though approximately 4000 m 3 /s of water flows through the Magdalena River in the study area, it is classified as a poor water quality source based on the Water Quality Index (WQI) 47 .
The drainage network discharges from the South-East to the North-East direction in some swamps or in the Magdalena River (Fig. 1b).Historically, water supply for human consumption in the rural area of VMM has been sourced from surface water.However, in recent decades, there has been a shift in the water usage patterns associated with a reduction of surface water availability and water contamination.In the study area, it has been identified that only the municipalities of Aguachica, San Martín, and Barrancabermeja are equipped with wastewater treatment plants for the treatment of both domestic and non-domestic wastewater, thereby facilitating compliance with Colombian regulations for discharges 48,49 There are no rural settlements with wastewater treatment systems.The disposal of domestic wastewater in rural settlements often occurs in public spaces, water bodies, or directly onto the soil, with excreta primarily managed through handmade pit latrines. www.nature.com/scientificreports/

Hydrogeological context
During the Mesozoic, the central and eastern Andes Mountain ranges were elevated, resulting in the formation of the MMV basin.It consists of sedimentary strata from the Cretaceous Period over a Jurassic igneous metamorphic basement with an average thickness of 8500 m [50][51][52][53] .To the south, there are outcrops of the Real Group, a Neogene formation whose thickness varies from 450 to 3500 m.In the north, there are Quaternary deposits outcrops underlying Real Group.Quaternary deposits have been categorized into eight geoforms, including alluvial fans and terrace deposits in the northeast, and others influenced by the deposition processes of the Magdalena River, Lebrija River, Sogamoso River, and primarily swamp zones, such as fluviolacustrine deposits, alluvial terraces, channel fluvial deposits, alluviums, terrace and deltaic cone deposits, and floodplain deposits 53,54 .At the regional scale, groundwater flow is governed by the topographic effects 55,56 , following in the South-North direction discharging in the swamps system of the Magdalena river.Moreover, there are secondary flows in the East-West direction, coming from the Eastern range and discharging into the Magdalena River.

Water sampling
88 surface water samples and 428 groundwater samples were collected in three (3) field sampling campaigns for the analysis of physicochemical and microbiological parameters.The first campaign was conducted during the dry period in February 2020 (Fig. 2a); the second was carried out in November 2020 during the wettest period of the hydrological (Fig. 2b) year and, the third one in March 2021 at the end of the dry period (Fig. 2c).2020 was cataloged as the La Niña phenomenon year, characterized by the increase in precipitation events.
The selection of groundwater sampling points was based on five criteria: (1) Spatial Distribution: Groundwater points were selected to ensure representativeness in the Neogene-Quaternary aquifer units outcropping in the study area, across different land uses/land covers, and within the main settlements of the study area.The selection of points less than 1 km apart was avoided to enhance spatial distribution.(2) Depth Distribution: Wells with depths greater than 10 m were chosen to ensure representation across the various aquifer units.(3) Accessibility: Selected points are within 1 km of the Colombian road network (primary, secondary, or tertiary) to minimize logistical issues and ensure sample integrity.However, this represents a limitation in the analysis of samples in the

Statistical analysis
Statistical analyses were performed using Python-Stats package.Basic statistics (mean, standard deviation, minimum and maximum values) were estimated for each parameter.The Kruskal-Wallis non-parametric statistical test was used to evaluate significant differences among samples, sampling period, and location 59,60 .This test enables the comparison of more than two independent samples, whether they have the same or different sample sizes, and mathematically, it is expressed as 59,61 : (1) where N corresponds to the total number of observations across all groups, g is the number of groups, n i is the number of observations in group i , r ij represents the rank of observation j from group i , is the average rank of all observations in group i , and r * = 1 2 (N + 1) is the average of all the r ij values 61 .

Water quality evaluation and spatial analyses
Water quality assessment followed the Colombian regulations 43 and US EPA standards for human consumption 42 as most people take water directly from wells with no treatment.Nevertheless, as Colombian regulation is less restrictive, the results shown here refer to US EPA standards.As a result, the portion of water samples above the allowed limits per parameter was quantified and identified along the study area.Their spatial behavior was compared with the land-use cover maps (1:100,000) 62 and with the location of the most crowded settlements.The Land Cover geographic object for the period 2018 provided by the Agustin Codazzi Geographic Institute (IGAC) at a scale of 1:100,000 62 was interpreted by means of Landsat 8 remote sensing data.Based on the experience and scope of the research, the 36 land cover classes proposed by IGAC for the study area were reclassified into 5 categories: (1) agriculture, (2) urban fabric, (3) water bodies, (4) pastures, and (5) forests.From this classification, the land use/land cover type was determined at each sampled point.

Human health considerations
The occurrence of total and fecal coliforms in water intended for human consumption has the potential to produce ADD, which is considered the leading cause of death in children under 4 years of age in developing countries 21 .For this reason, this research analyzed under-four children mortality caused by ADD during the year 2019, based on the information reported by the National Institute of Health of Colombia through the Public Health Surveillance System (SIVIGILA) for municipalities located in the study area 63 .Municipalities with mortality in children under four years of age were identified and related to the presence of pathogens in water supply points within their jurisdiction.
In the MMV, some settlements take groundwater for their water supply systems but lack drinking water treatment.The same occurs in the farms, where groundwater supplies agricultural activities and human consumption.
As the healthcare system in the rural area of VMM is precarious or non-existent, cases of ADD are often referred to medical centers in nearby major cities, some of which are even located outside the study area.Centralized infant mortality data are based on reports from clinics and hospitals, not on municipal reports, which could introduce biases in the information.Therefore, a risk analysis of water consumption in the Middle Magdalena Valley was conducted based on the laboratory analyses.

Colombian Drinking-Water Quality Risk Index (CDWQRI-IRCA)
One method for assessing the risk related to the quality of water designated for human consumption involves the calculation of the Water Quality Indices (WQI) [64][65][66] .These indices enhance the understanding of water quality and serve as a tool for decision-making in the absence of specialized knowledge in specific fields related to water resources 67 .In this research, the Colombian Drinking-Water Quality Risk Index (CDWQRI), known in Colombia as IRCA, was estimated.This index is an arithmetic measure that relates the sum of weighted values, known as risk scores, for each of the parameters considered in the analysis, which do not meet the standards proposed by Colombian regulations 43 , to the total weighted values of all measured parameters, including those that meet the regulations.The CDWQRI-IRCA for each analyzed sample is determined by the Eq.(2)   where RSUAP represents the risk score assigned to unacceptable parameters, while RSMP corresponds to the risk score of measured parameters.The weighted values are presented in the Table 1.
The risk level for human health evaluated using the CDWQRI-IRCA is categorized into five classes: Unhealthy, High, Medium, Low, and without risk.Higher values are associated with water unfit for human consumption, while values close to zero are considered suitable for human consumption Table 2.

Water quality in the MMV
The Kruskal and Wallis non-parametric test 59,60 revealed statistical differences among the sampling period (wet or dry) and the location of the water sample (referred to the northern and the southern area) with p-values below 0.05 for Ca, Mg, HCO − 3 and Cl − .Accordingly, most results are presented by sampling period.In Table 3, mean, minimum, maximum, and standard deviation values associated with both the physicochemical and microbiological analyses are presented for each sampling period, and Table 4 displays an overall statistical analysis for all the analyzed samples (Supplementary Information).
The percentage of water samples that either comply or do not comply with the US EPA standards 68,69 and Colombian regulation 43,49 for human consumption is depicted in the Fig. 3.Among the 28 parameters evaluated, the Colombian regulations are found to be more stringent than the EPA standards for 10 specific parameters (Ba, Cyanide, Zn, Cu, Cr, F − , Hg, NO − 2 , Pb, and Se).Conversely, for three parameters (Mn, Cd, and pH), the Colombian standards are less stringent.Fecal coliforms and Total coliforms exceed the quality limits in 58% and 89% of the analyzed samples, respectively, in relation to both regulations.pH, color, Fe, Mn, NO Overall, total coliform occurrence is observed throughout the study area in the majority of the sampled sites (at least 89% of the analyzed samples), with higher concentrations near major population centers, except in the vicinity of Barrancabermeja, a major city with a robust sewer system in its urban area (Fig. 4).The most significant concentrations of pollution attributable to the presence of total coliforms were observed in the surface water samples under analysis, a finding consistent with expectations given that both domestic and industrial wastewater discharges are typically released into water bodies without prior treatment.Among the various surface water sites monitored, the Quebrada Cayumbita, situated in the southern region of the Sabana de Torres municipality, exhibited the highest level of total coliform contamination, registering 420,000 MPN/100 mL.In contrast, groundwater samples revealed the greatest total coliform concentration in a well located within the urban perimeter of the Aguachica municipality.Despite possessing a drinking water treatment system that is supplied from a surface water source, Aguachica faces water supply issues, compelling the local populace to resort to groundwater as a supplementary resource 70 .
Analyzing histograms with smoothed function fitting for total coliform analyses in each field campaign, a higher number of samples exhibiting total coliform presence was observed during the wet period (Fig. 4).The behavior associated with the occurrence of total coliforms in the analyzed samples is consistent across all three campaigns.The histograms indicate that the occurrence of total coliforms in the analyzed samples is concentrated at values below 50,000 MPN/100 mL, consistent with the analysis of central tendency measures, which reveal a mean of 13,136 MPN/100 mL and a standard deviation of 41,450 MPN/100 mL.
During the fieldwork, it was observed that the sampled wells, which are mostly located in scattered rural households, are situated near handmade pit latrines without proper hygienic protection or septic tanks and, in Table 1.CDWQRI-IRCA risk score.Risk scores associated with each physicochemical and microbiological parameter used to estimate the CDWQRI-IRCA are presented 43 .www.nature.com/scientificreports/areas where backyard livestock production takes place.Consequently, contamination of groundwater sources by both total and fecal coliforms is expected, as evidenced by the laboratory results obtained.
Nitrates exceeded the permissible limit (10 mg/L) in more than 10% of the water samples.The spatial distribution of nitrates shown in Fig. 5 revealed higher concentrations in urban areas of Aguachica, Sabana de Torres and San Martin, probably as a result of nitrogen-rich effluent discharges counting for the lack of proper sewage treatment.High values were observed too, in specific agroindustrial zones near Puerto Wilches, where there are extensive areas of African oil palm crops, and towards Barrancabermeja, an area that was previously dedicated to crops such as rice and cotton but that nowadays is being converted into livestock [71][72][73] .
Histograms with smoothed function fitting for nitrate analyses indicate that nitrate occurrence is unrelated to the hydrological period during sampling Fig. 5.This is further substantiated by analyzing measures of central tendency, revealing consistent mean values of 4.3 mg/L for both dry and wet periods, with standard deviations of 9.4 mg/L and 10.7 mg/L, respectively.This finding is unexpected, considering that in Colombia, traditional farming practices often involve applying larger quantities of fertilizer directly to the soil, primarily during wet periods, to mitigate the effects of fertilizer runoff.
Additionally, we conducted a spatial analysis for pH, iron and manganese, parameters that did not comply with the US EPA standards (Fig. 6).pH exceeded the permissible limit in 195 samples (40%).These samples are primarily situated between the Lebrija and Sogamoso rivers, an area predominantly characterized by African oil palm crops (Fig. 7) Iron was exceeded in 65 water samples (16.7%), mainly in the southern of the study area, possibly associated with the presence of soils rich in iron oxides.Regarding manganese, a noteworthy 21.4%, equivalent to 80 out of the total water samples were above the permissible limit.Elevated manganese concentrations are widely dispersed throughout the entire study area, implying a connection to inherent natural factors tied to the mineralogical traits of exposed geological formations, rather than anthropogenic influences 74,75 .
Finally, the presence of metals such as Ag, Pb, As, Ag and Hg is an alert for human consumption, even though those are reported in less than 3% of the water samples.Table 3. Descriptive statistics on the 28 physicochemical and microbiological parameters of the water quality analyzed in the MMV.These are presented according to the hydrological period in which the sampling was conducted.

Effect of land-use/landcover on groundwater quality
Land cover in the study area was classified into six (6) categories grouping different land use/landcover characteristics.The analysis revealed that 43.4% of the surface is covered by Grasslands, which is the predominant category and is used for extensive livestock farming.Agriculture lands occupy the second largest area with 24.9%, with African palm plantations for agro-industry and crops such as rice and corn being the most common.Forests account for 14.7% of the surface.Wetlands cover 7.9% of the area, including all marshy areas and even zones with aquatic vegetation on water bodies.Water bodies, including both lentic and lotic surface water bodies, cover 7.5% of the study area.Urban areas represent the smallest proportion of land use/land cover in the MMV, estimated to be less than 1.7% of the total surface area.Impacts of land use and land cover on groundwater quality vary depending on the specific land categories under consideration.Agricultural lands are mainly located in the southern portion of the study area, spanning between the municipalities of Puerto Wilches and Barrancabermeja.In the northern part of the study area, grasslands predominate, primarily used for livestock farming, interspersed with scattered agricultural plots throughout the region.The precipitation distribution could potentially be linked to this behavior, given that the southern sector receives a total annual rainfall of 4000 mm, while in the northern falls below 1000 mm.
Figure 8 illustrates the contrast between land use/land cover and the spatial distribution of water samples analyzed by hydrological period.The samples are categorized into three groups: (1) samples in which all the analyzed parameters are below US EPA limits, (2) those samples in which one or two parameters fail to meet the standards, and (3) those samples in which three or more parameters exceed the US EPA standards.It was observed that during the dry periods, the analyzed sampled sites complied with a greater number of parameters, especially in the northern area, characterized by grasslands.On the other hand, in the southern area, a greater number of sampled sites exhibited more than three parameters exceeding the US EPA standards in both during dry and wet periods.
Regarding the spatial distribution of nitrates in groundwater, higher concentrations were found near the main population centers lacking proper sewage systems (Aguachica, San Martín, and San Alberto), suggesting a potential association with wastewater discharges (Fig. 5).On the other hand, in the central area of the study area, characterized by scattered rural dwellings and small, low-density populations, the lowest nitrate concentrations were detected (< 5 mg/L).Temporally, the highest nitrate concentration rates occurred during the wet period, which could be linked to the increased use of pesticides, typical in Colombian agriculture practices 76 .Lastly, shallow groundwater acidification was observed particularly concentrated between the Sogamoso River and the Lebrija River (Fig. 7).This area is characterized by African oil palm cultivation, suggesting a potential relationship between these crops, soil acidification, and consequently, groundwater acidification 77,78 .

Assessment of human health risks
The risk of disease occurrence or health-related alterations associated with water consumption was assessed using the CDWQRI-IRCA.This index is calculated by considering the non-compliance with physicochemical and microbiological water parameters in comparison with the Colombian regulation for drinking water.
The CDWQRI-IRCA was computed for 458 water samples, spanning values from 0.0 to 79.2, with an average of 50.9, a median of 61.5, and a standard deviation of 21.1.Among the 458 samples, 386 of them, accounting for 84.28% of the total analyzed, were classified as high risk, 4.59% indicated a medium risk, 5.46% a low risk, whereas 5.68% exhibited no risk for human consumption.It is crucial to highlight that none of the samples were deemed unsuitable for consumption based on sanitary considerations.
From CDWQRI-IRCA histograms, no significant temporal trend was observed that would allow correlating the hydrological period with the risk to human health due to the consumption of groundwater (Fig. 9).During the dry periods of 2020 and 2021, it was estimated that 90.87% and 71.84% of the sampled water points were categorized as high risk, respectively, and during the rainy period of 2020, this figure was 86.16% (Table 5).
Spatially, in the rural area characterized by dispersed rural housing and therefore by groundwater points destined for domestic use for individual housing units, a noticeable trend was identified in the samples characterized with a high-risk level (Fig. 9).In the vicinity of the main populated centers such as Barrancabermeja, San Martín, Sabana de Torres, Aguachica, San Alberto, and San Martín, sampled sites did not exhibit a high level of risk, including a few associated with municipal water supply systems.
Therefore, a detailed analysis of the points associated with water supply systems, both municipal and rural, was conducted.Sampling took place at storage tanks to estimate the actual water quality consumed by the population.It was found that 57 points, which account for 76%, were categorized with a high-risk level for human consumption.10 points (13%) are at a medium or low-risk level, while only 8 points (11%) pose no risks to human health.

Direct impacts on human health
The National Institute of Health of Colombia reports two (2) deaths among children under 4 years due to ADD in 2019.One in Cantagallo and another in San Pablo on the western bank of the Magdalena River which may www.nature.com/scientificreports/be related to the fact that 84% of the analyzed samples were classified as high-risk for human consumption.Unfortunately, because of the absence of health centers in rural areas and the lack of a systematic record of ADD cases, it is not possible to correlate the consumption of polluted water with the occurrence of diseases or hazards to health.In fact, it seems that people have created resistance to the DAA as not native people used to get sick when visiting the region and consuming raw water.

Conclusions
This study assessed groundwater quality in the Middle Magdalena Valley in Colombia through the physicochemical and microbiological analysis of 28 water quality parameters.The research was conducted in three campaigns at 190 groundwater sampling points.The results showed that Ag, Se, Hg, F − , Zn, and Ba met the water quality standards proposed by the US EPA for human consumption in all analyses.Among the other parameters analyzed, both fecal and total coliforms exceeded the quality standard in a higher percentage of samples, specifically 58% and 89%, respectively.Given that the majority of the sampled groundwater points were located in rural dwellings, where groundwater is typically used for human consumption without any treatment, and this behavior is widespread in the rural area of the Middle Magdalena Valley, an assessment of human health risks associated with the consumption of these waters was conducted.This assessment was carried out using the CDWQRI-IRCA.
The results obtained from the CDWQRI-IRCA evaluation indicate that only 5.7% of the analyzed samples pose no health risk, while 84.3% of them represent a high risk, demonstrating a correlation with the results of the total coliform analysis.As a consequence of the poor quality of groundwater, often consumed by the rural population of the MMV, two deaths of children under 4 years old due to Acute Diarrheal Disease were recorded in 2019.The study identified that the main problem concerning water quality and its public health implications is biological contamination due to the occurrence of fecal and total coliforms at the measured groundwater points.Based on these findings, it is recommended to implement drinking water treatment systems, with an emphasis on pathogen removal and prioritizing their establishment in rural areas.Additionally, it is recommended to assess the feasibility of implementing low-cost household systems for water purification, including community aqueducts, in dispersed rural communities, with the support of organized grassroots communities.These recommendations align with the 2030 Agenda and the Sustainable Development Goals-An opportunity for Latin America and the Caribbean 79 .However, this poses a significant challenge in rural Colombia, where, by the end of 2023, only 57.8% coverage has been achieved in populated centers and scattered rural areas 80 .
While the previously mentioned solution represents a structural and largely definitive answer to the problem, alternative measures are needed to address this issue.Among these, the need to ensure protection systems in wells to be constructed, such as sanitary seals, stands out; avoiding the construction of wells in the vicinity of handmade pit latrines or areas of wastewater discharge, especially if they are upstream of the wells to be built for human supply; abandoning wells that are within 50 m of wastewater discharge points or handmade pit latrines 81 ;   intensifying governmental efforts in constructing effective wastewater treatment systems, both in urban and rural areas; ensure compliance with national regulations for the discharge of industrial wastewater by constructing industrial wastewater treatment systems in agro-industrial companies within the territory; establish systematic groundwater monitoring plans at both national and regional levels, incorporating participative monitoring processes and promote water governance; and implement environmental education campaigns and processes designed to enhance awareness and understanding of water conservation, utilization, and management.These initiatives should be aligned with a National Environmental Education Policy and thus be integrated into both national and regional environmental education programs to ensure the sustainability of the actions undertaken.
This research represents an effort to assess the quality of groundwater and understand its potential public health impacts.However, there is a need for greater understanding of the possible impacts that climate change could induce, particularly regarding changes in land use, anthropogenic pressures on groundwater, and aquifer recharge, which are intrinsically related to alterations in precipitation patterns.Therefore, it is proposed that long-term groundwater quality monitoring strategies be implemented and that future studies be conducted to examine the effects of climate change on groundwater quality in the MMV at a regional scale.

Figure 1 .
Figure 1.Study area showing (a) general location; (b) elevation in meters above mean sea level, municipalities, and major surfaces water bodies.This figure illustrates the spatial location of the MMV, with a particular focus on its topography.The HydroSHEDS DEM (World Wildlife Fund, 2006-2013), generated by NASA's SRTM (Shuttle Radar Topography Mission) at a spatial resolution of 92 m, was employed for this purpose; (c) spatial distribution of mean annual precipitation.Additionally, this figure depicts the temporal distribution of precipitation evaluated at three meteorological stations managed by the Institute of Hydrology, Meteorology, and Environmental Studies of Colombia (IDEAM).Figures were generated using ArcGIS Pro software, Version 3.1 (https:// pro.arcgis.com/ es/ pro-app/3.1/ get-start ed/ whats-new-in-arcgis-pro. htm).

Figure 2 .
Figure 2. Map showing groundwater sampling sites during the (a) dry period in March 2020; (b) rainy period in November 2020; and (c) dry period in March 2021.These are categorized by water source type, including springs, shallow wells, wells, rivers, and swamps.

Figure 3 .
Figure 3.Comparison of water quality analyses between the standards proposed by the US-EPA and the Colombian Regulations for Drinking Water.Green bars represent the percentage of analyzed samples that meet the quality standards for human consumption.Orange bars indicate the percentage of samples that fail to meet the standards.Blue bars correspond to the percentage of samples from which no conclusions can be drawn.

Figure 4 .
Figure 4. Spatial distribution of total coliform occurrence during (a) the dry period of 2020; (b) the rainy period of 2021; and (c) the dry period of 2021.Additionally, histograms with smoothed function fitting for the total coliform analyses evaluated in each field campaign are presented in (d).

Figure 5 .
Figure 5. Spatial distribution of nitrate occurrence evaluated from the three conducted sampling campaigns.(a) Corresponds to the dry period of 2020; (b) to the rainy period of 2021; and (c) to the dry period of 2021.In (d), histograms with a smoothed function fit are presented for each field campaign.

Figure 6 .
Figure 6.Spatial distribution of water points where analyzed samples surpass the US EPA thresholds.In (a) the analysis is presented for pH; in (b) for iron; and in (c) for manganese.

Figure 7 .
Figure 7. Spatial distribution of pH based on the analyzed samples.(a) Presents the analysis for the dry period (2020); (b) for the rainy period (2021); and (c) for the dry period (2021).pH values fluctuated between 3.57 and 8.83 in the analyzed samples.

Figure 8 .
Figure 8. Land cover/soil use map in the MMV and the number of parameters above the US EPA standard per sample point.Land cover/use has been categorized into urban areas, wetlands, surface waters, grasslands, crops, and forests.Additionally, out of the 28 physicochemical parameters analyzed, the quantity of parameters exceeding the thresholds established by the US EPA for each sampled water point is presented.This analysis was conducted for the three campaigns: (a) dry period (2020); (b) rainy period (2021); and (c) dry period (2021).

Figure 9 .
Figure 9. CDWQRI-IRCA values for the analyzed hydrological periods.The results of the CDWQRI-IRCA have been categorized as unhealthy, high risk, medium risk, low risk, and no risk for each conducted sampling campaign: (a) dry period (2020); (b) rainy period (2021); and (c) dry period (2021).Additionally, in (d), histograms with smoothed function fitting are presented for the results of the CDWQRI-IRCA analyses, discriminated by sampling campaign.

Table 4 .
Descriptive statistics of the physicochemical and microbiological analyses without discriminating the hydrological sampling period.

Table 5 .
Assessment of human health risk based on CDWQRI-IRCA for the analyzed hydrological periods.